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P-t[ ABSTRACT 

q , We present high spatial resolution optical imaging and polarization observations of the 

?H ■ PSR B0540-69.3 and its highly dynamical pulsar wind nebula (PWN) performed with Hub- 

c/3 \ ble Space Telescope, and compare them with X-ray data obtained with the Chandra X-ray 

d . Observatory. In particular, we have studied the bright region southwest of the pulsar where 

a bright "blob" is seen in 1999. In a recent paper by De Luca et al. it was argued that the 
| "blob" moves away from the pulsar at high speed. We show that it may instead be a result of 

local energy deposition around 1999, and that the emission from this then faded away rather 
than moved outward. Polarization data from 2007 show that the polarization properties show 
, dramatic spatial variations at the 1999 blob position arguing for a local process. Several other 

r**"* ■ positions along the pulsar- "blob" orientation show similar changes in polarization, indicat- 

ed ' ing previous recent local energy depositions. In X-rays, the spectrum steepens away from the 

qq . "blob" position, faster orthogonal to the pulsar- "blob" direction than along this axis of orien- 

■ tation. This could indicate that the pulsar- "blob" orientation is an axis along where energy in 

| the PWN is mainly injected, and that this is then mediated to the filaments in the PWN by 

7— I • shocks. We highlight this by constructing an [S II]-to-[0 III]-ratio map, and comparing this 

to optical continuum and X-ray emission maps. We argue, through modeling, that the high 
[S II]/[0 III] ratio is not due to time-dependent photoionization caused by possible rapid X- 
ray emission variations in the "blob" region. We have also created a multiwavelength energy 
spectrum for the "blob" position showing that one can, to within 2<r, connect the optical and 
X-ray emission by a single power law. The slope of that power-law (defined from F v = v~ c ' v ) 
would be a Y = 0.74 + 0.03, which is marginally different from the X-ray spectral slope alone 
with a v - 0.65 + 0.03. A single power-law for most of the PWN is, however, not be possi- 
ble. We obtain best power-law fits for the X-ray spectrum if we include "extra" oxygen, in 
addition to the oxygen column density in the interstellar gas of the Large Magellanic Cloud 
and the Milky Way. This oxygen is most naturally explained by the oxygen-rich ejecta of the 
supernova remnant. The oxygen needed likely places the progenitor mass in the 20 - 25 M 
range, i.e., in the upper mass range for progenitors of Type IIP supernovae. 

Key words: pulsars: individual: PSR B0540-69.3 - ISM: supernova remnants - ISM: SNR 
0540-69.3 - supernovae: general - Magellanic Clouds 
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1 INTRODUCTION 

PSR B0540-69.3 is a 50.2 millisecond pulsar in the Large Mag- 
ellanic Cloud (LMC). A detailed review about the discovery 
and observations of the pulsar, its wind nebula, as well as the 

* E-mail:natalia@astro. su.se 



whole supernova rem nant in all wavelength bands is provided by 
IWiriiamsetaT1 ( l2008h . 

PSR B0540-69.3 is often referred to as the "Crab twin" be- 
cause it is of similar age (~ 1000 years), it has a similar pulse pe- 
riod and that it has a pulsar wind nebula (PWN) surrounding it, sim- 
ilar to that of the Crab. It has also been suggested that the detailed 
structures of the two PWNe, as revealed for PSR B0540-69.3 in 
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X-rays dGotfhelf & WansfcOOd) , are similar, and that PSR B0540- 
69.3 may have a jet and a torus like the Crab. For further com- 
parison between the two pulsars and their PWNe we refer to 
ISerafimovich et al. I Jiooj) . 

There are, however, important differences between the two 
pulsars and their surrounding nebulae. The most obvious is that 
the supernova ejecta of SNR 540-69.3 are oxygen-rich, although 
with some hydrogen mixed in dSerafimovich etal. H 2005). making 
it likely that SNR 0540-69.3 and PSR B0540-69.3 (together hence- 
forth referred to as "0540") stem from a Type IIP sup ernova ex- 
plosion of a massive (~ 20 Mo) star dChevaliedl2006l) . whereas 
the Crab pro genitor was a much less massive (8-10 M Q ) star (e.g., 
lHestej2008h . In addition, 0540 has a normal supernova shell of fast 
ejecta, something which h as not yet been fully confirmed to exist 
around the Crab (see, e.g.. lTziamtzis et al. 2009). This could mean 
that the Crab ejecta have much less kinetic energy than the ejecta 
of 0540 which appears to be normal in that respect. 

Closer to the center, 0540 appears much more asymmetric 
than the Crab Nebula, with much of the emission c oming from 
a regi on a few arcseconds southwest of the pulsar. iMorse etal. I 
(2006) show that this asymmetric appearance, as well as red- 
ward asymmetry in i ntegrated profiles of emis sion lines (see also 
llCirshner et al. | [l989: Serafi movich etal .1120051) , does not signal an 
overall asymmetry of the PWN as the maximum velocity of emis- 
sion lines toward and away from us are nearly equal. 

Even closer to the center, a detailed comparison between the 
PWNe of the two pulsars is more di fficult due to th e large dis- 
tance to PSR B0540-69.3 (~ 51 kpc. |Panagia|[2 005) as opposed 
to the 2 ± 0.5 kpc to the Crab dKaplan et al. 1120081) and references 
there in); the semi-major axis of the torus in the Crab would only 
subtend ~ V.'S at the distance of PSR B0540-69.3, and details 
such as the changing structure of the Crab wisps would have been 
close to impossible to disentangle even with the HST or ground- 
based adaptive optics. It therefore comes as no surprise that there 
is yet no counterpart for PSR B0540-69.3 and its PWN to the de- 
tailed synchronized optical and X-ray study of the Crab PWN by 
Hest er etal. ] J2002I) . These authors revealed the presence of out- 
ward moving equatorial wisps with velocities of ~ 0.5 the speed of 
light. Many time- varia ble subarcsecond s tructures e xist in the Crab 
PWN in the opti cal dHester et al. I [20021). near-IR dlvlelatos et al. I 
l2005t) and X-rays (Hester etal. ll2002UWeisskopf etal. Il2000l) . This 
is discussed in detail by Hestei 1 20081) . 



The first indication that also 0540 displays chan ges in its struc- 
ture ov er a relatively short time period was shown bv lDe Luca et al. I 
d2007l) . They found a 10 year flux variation (between 1995 and 
2005) in the southwest direction where the PWN has its strongest 
emission and suggested that they had detected a hot spot moving 
at ~ 0.04c. De Luca et al. noted that the hot spot could be simi- 
lar to a time-varying arc-like feature in the outer Crab Nebula, and 
that a pulsar jet in 0540 could be directed toward the bright south- 
w est region rather than per pendicular to this, as suggested earlier 
bv lGotthelf & Wand d200d) . 

In this paper we prese nt a study of th e PWN of PSR B0540- 
69.3 using the same data as lDe Luca et al. ( 2007) , complemented 
with all other HST data available for 0540 as well. In particular, 
we have included recent HST polarization data from 2007. We also 
include all available X-ray data of 0540. 

The paper is organized as follows: in Sect. 2 we describe and 
discuss the observations, in Sect. 3 the analysis and results, and in 
Sect. 4 we conclude with a discussion. 



2 OBSERVATIONS 

2.1 HST observations and data analysis 

The pulsar field was obseved with HST/WFPC2 in 2005 (Program 
ID 10601) using the two broadband filters F547M and F555W 
(see Table [TJ. As opposed to earlier HST/WFPC2 observations of 
PSR B0540-69.3 we used a dithering procedure for all images to 
better avoid cosmic ray contamination. PSR B0540-69.3 and its 
PWN were exposed on the PC chip, and the observations were 
made in a two-gyro mode. The CCD amplification, i.e., gain was 
7.12 and readout noise=5.24. The reason for choosing the contin- 
uum filter F547M was that the main goals of the observations were 
to study the proper motion of the pulsar, as well as the spectrum of 
the pulsar and the PWN. In earlier observations the pulsar and PWN 
were mainly studied in filters largely contaminated with line emis- 
sion from filaments in 0540. Continuum filters like F547M avoid 
this contamination. The 2005 data have been discussed by De Luca 
et al. (2007), and an upper limit to the proper motion of the pulsar 
was estimated, as well as identifying the apparent displacement of 
a feature in the PWN. We discuss this further here, as well as use 
the data for spectral studies of the PWN. The 

To study time variability of the PWN we also used available 
archival HST/WFPC and WFPC2 data obtained from 1992 to 2007 
(Table [Q. There are three different storage sites of HST observa- 
tions, namely the MASlQ, CADC0, and ST-ECr§ All three archives 
use slightly different calibration plans, or pipe-line reductions of 
raw data. We found that the pipe-line reduced data of the same 
WFPC2 observations from these archives give different results for 
the photometry of PSR B0540-69.3 and its neighborhood. In previ- 
ous work (Serafimovich et. al. 2004) we did not resort to pipe-line 
reductions, but made careful manual reductions of the raw data. Us- 
ing photometric results from that work we found that only the data 
from the MAST archive with their recent pipe-line calibration plan 
can reproduce our previous photometry. To avoid systematic differ- 
ences due to pipe-line effects in our analysis we therefore used the 
data only from the MAST archive with the most recent calibration 
plan for all epochs of HST/WFPC2 observations listed in TableQ] 

All data were reduced using I RAF in combination with STS- 
DAS and DITHEF0 packages. Pre-calibrated data include refer- 
ences to the latest distortion correction files available at the STScjJ 
Using STScI recommendation^ we used these files together with 
the Multidrizzle package for each individual image independently. 
Cosmic ray cleaning and combining of the images were excluded 
at this step. This exclusion was essential for the 1995 and 1999 
epochs. These observations were performed without dithering in 
between exposures, and after cleaning and combining, Multidriz- 
zle introduced extra noise in the pixel distribution which led to an 
overestimate of faint source fluxes. The overestimate for F547M 
in the 1999 epoch was found to be a factor of around two. After 
Multidrizzle, the individual images were converted from the unit 
of electrons per second to counts using the IRAF/imgtools/imcalc 
utility, e.g., (image * exposure)/ gain. Final image combining was 
performed using the standard IRAF imcombine utility with cosmic 
ray rejection by setting reject=crreject. 

The data obtained in 2005, and later, clearly show that 
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Table 1. HST/WFPC2 optical observations of the PSR B0540-69.3 field used in our study. 



Date 


Total 




Filter 




Proposal 


Comments 




exp." 


name 


PHOTPLAM* 


PHOTBW r 


ID 




17Nov92" 


2x400 


F547M 


5483.9 


205.5 


4244 


continuum 


190ct95 


2x300 


F555W 


5442.9 


522.2 


6120 


overlaps with [O III] A 5007 


170ct99 


2x300 


F336W 


3359.5 


204.5 


7340 


continuum 


170ct99 


8x1300 


F502N 


5013.3 


48.3 




centered on [O III] A 5007 


170ct99 


2x400 


F547M 


5483.9 


205.5 




continuum 


170ct99 


6x1300 


F673N 


6732.3 


30.7 




centered on [S II] AA 6716,6731 


170ct99 


2x200 


F791W 


7872.5 


519.8 




continuum 


15Nov05 f ' 


4x260 


F547M 


5483.9 


205.5 


10601 


continuum 


15Nov05 f ' 


3x160 


F555W 


5442.9 


522.2 




overlaps with [O III] A 5007 


21Jun()7 


3x260 


F336W 


3359.5 


204.5 


10900 


continuum 


21Jun07 


3x260 


F450W 


4557.3 


404.2 




continuum 


21Jun07 


3x100 


F555W 


5442.9 


522.2 




overlaps with [O III] A 5007 


21Jun07 


3x140 


F675W 


6717.7 


368.3 




overlaps with [S II] AA 6716,6731 


21Jun()7 


3x200 


F814W 


7995.9 


646.1 




continuum 



" Total exposure is given in number of images times the * Pivot wavelength of the filter band measured in A. 

exposure time of individual images in seconds. c Width of the filter band measured in A. 

d The 1992 observations were obtained with HST/WFPC e Our observations. 

instead of HST/WFPC2. 
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Table 2. Chandra X-ray observations of the PSR B0540-69.3 field used in 
our study. 



Ill 1U T1X Tlj? liX IIMU 

at 



Figure 1. A 15"xl5" image of the field around PSR B0540-69.3 obtained 
in the F547M band with HST/WFPC2 in 2005 (Table [T}. The pulsar was 
exposed on the PC chip and its position is marked by an arrow. The diffuse 
emission surrounding the pulsar is the pulsar wind nebula (PWN). A highly 
variable area or "blob" is marked by a circle and an arrow. Two typical 
traces of the CTE effect are marked by arrows (see text for more details). 



the Charge Transfer Efficiency (CTE) effect dRiessI feOOCh l on 
HST/WFPC2 images was getting worse with time (see Figure [TJ. 
The CTE effect is very complicated to compensate for as it is time, 
position and flux dependent, i.e., it is unique for every particular 
point source and almost unknown for extended objects. In order to 
more reliably study the structure of the PWN, we used the multi- 
scale filtering method described in Sect. l2.5l and l3.ll 



Date 


Exposure 
(ks) 


Instrument 


ObsID 


31Aug99 


17.79 


HRC-I 


132 


21Jun00 


10.04 


HRC-I 


1735 


21JunOO 


10.04 


HRC-I 


1736 


22Jun00 


10.07 


HRC-I 


1741 


15Feb06 


40.36 


ACIS-S 


5549 


16Feb06 


39.98 


ACIS-S 


7270 


18Feb06 


38.56 


ACIS-S 


7271 



observations of 0540 were obtained as part of the instrument cali- 
bration plan in 1999 and 2000, only a few months before and after 
the HST 1999 observations described in TableQ] They contain very 
limited spectral information, but provide the best spatial resolution 
(pixel size (J.' 132) in X-rays, and are therefore the most useful for 
the study of the PWN morphology. The ACIS-S observations (PI 
S. Park) have lower spatial resolution (pixel size C.'492), but pro- 
vide better spectral information. The ACIS-S observations started 
only three months after our HST observations in 2005, and are use- 
ful for an optical/X-ray comparison, assuming the PWN structure 
did not change substantially in three months. Another advantage of 
these observations is that the 1/4 ACIS-S subarray mode was usecQ 
and the data are much less contaminated by the CCD pileup effects 
from the high count rate of the pulsar in the center of the PWN, 
as opposed to pr evious ACIS obser vations, e.g., the data obtained 
in Nov 23 1999 jPetre et al. I J2007t) . Obs ID 119, exposure 28.16 
ks). For the latter, significant pileup precludes a reliable spatial and 
spectral analysis of the central PWN regions. The HRC-I data are 
not affected by the pileup. 



2.2 Chandra X-ray observations 

To compare the structure of the PWN in the optical and X-rays, 

we utilized the Chandra archival data listed in Table [2] The HRC-I 7 The frame lime of the 1/4 mode is < 1.1 s. For details see asc.hiirviird.edii/proposer/POG/himl/^hapfi.hlml 
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Table 3. HST/WFPC2 polarization observations of the PSR B0540-69.3 
field used in our study. 



Date 


Total 


Mode of obs. 


Pol. 


Proposal 




exp." 




angle 


ID 


05Nov()7 


600x3 


F606W+POLQ 


0° 


10900 


25Sep07 


600x3 


F606W+POLQ 


45° 


10900 


26Jul07 


600x3 


F606W+POLQ 


90° 


10900 


21Jun07 


600x3 


F606W+POLQ 


135° 


10900 



" The total exposure is given in number of images times 
the exposure time of individual images in seconds. 



2.3 HST/WFPC2 polarization observations 

To better understand the highly variable structure in the southwest- 
ern part of the 0540 PWN we utilized HST/WFPC2 archival polar- 
ization data. 

PSR B0540-69.3 and its wind nebula were observed in polar- 
ized light with HST/WFPC2 in 2007 (Program ID 10900) using 
the broad band filter F606W in combination with the POLQ po- 
larizer (see Table [3} . The field was exposed only on the PC chip 
during all sets of observations using different rotations of the HST 
to reproduce the required polarization angle (see column 4 in Ta- 
ble[3j. This gives the presently highest spatial resolution available 
for polarization observations, i.e., 070455 per pixel. Unfortunately, 
the pulsar and its nebula were placed mainly in the vignetted area 
of the PC chip. Since the part of the nebula was never more than 
4" into the bad area, the loss of light was not so big. However, this 
made calibrations slightly more cumbersome. 

Due to the long exposures (600 s) each frame exhibits thou- 
sands of cosmic ray (CR) hits. Most of these we removed by tem- 
poral median filtering of the three images, and the remaining hits 
are removed by la-cosmic which is included in the IDL software. 
The images were taken at the roll angles 0, 45, 90 and 135 degrees, 
and to align the images three of them must be rotated and shifted. 
The de-rotation of 45 and 135 degrees requires interpolation, and so 
do the shifts. We used a cubic spline interpolation scheme and even 
though we succeed in overlapping the images to an accuracy of 10 
milliarcseconds, it is clear that the resulting PSF (point spread func- 
tion) for point sources may be affected. This means that we must be 
careful in our interpretation of polarization close to point sources. 
On the other hand, the extended emission, which is the focus of the 
present investigation, is not affected by the interpolations. 

In this set of observations, the main source of instrumental 
polarization is the WFPC2 pick-off mirror. We have followed the 
calibrati on scheme as outlined in WFPC2 Polarization Calibration 
Theory dBiretta & McMasteill 19971) to calculate the Mueller matrix 
for 605 nm, and as a result we get following relations 

1 = 1.11* 

Q= 1.1 (Q obs + 0.065 l obs ) 
U= 1.6 (LU S + 0.034 I o/ „), 

where I, Q and U are standard Stokes parameters, and \ obs , Q oA , 
and U„/„ are Stokes parameters before correction for instrumental 
polarization. Test measurements of stars in the field show a small 
tendency of a remaining systematic polarization of ~ 2%, which is 
a bit more than expected. According to a detailed study by Jeffery 
(1991) of SN 1987A, which is located quite close to PSR B0540- 
69.3, the foreground polarization is only ~ 0.4%, so the ~ 2% sys- 
temic uncertainty we find probably also represents the uncertainty 
of our results for 0540. This is small in comparison to the observed 



polarization which is typically 10 - 30%, and therefore does not 
affect our results noticeably. 



2.4 Astrometric referencing 

Astrometric referencing of the HST/WFPC2 images was performed 
based on the positions of the astrometric standards selected from 
the USNO-B1 astrometric catalog A dozen of the USNO-B1 ref- 
erence objects can be identified within the PC2 chip FOV. Their 
lcr coordinate uncertainties vary between 50 and 900 mas with an 
rms value of about 330 mas. Pixel coordinates of the standards 
were derived using of the IRAF task imcenter with an accuracy 
of 0.05-0.1 PC2 pixel size (07042). The IRAF tasks ccmap/cctran 
were applied for the astrometric transformation of the images. We 
consequently discarded standards which are either over saturated or 
have relatively large catalog and image positional uncertainties. As 
a result, only seven stars were selected for the final astrometric fit. 
Formal rms uncertainties of the fit were best for the F814W filter 
ARA<0732 and ADec<0743, and the fit residuals are <076, which 
is compatible with the maximum catalog position uncertainty of 
the finally selected standards (0733). The rest PC2 images were 
referenced to the F814W image with the accuracy better than 0.1 
(07004) of the PC2 pixel size using several unsaturated stars from 
the PWN neighbourhood. Finally, the first F547M images obtained 
in 1992 with the WFPC at a worse spatial resolution (pixel size 
07 1) were aligned to the WFPC2 images with the accuracy of about 
0702. 

Inspection of the Chandra/HRC and ACIS images obtained 
in 2000 and 2006 showed that within a nominal Chandra pointing 
accuracy of <0.5 arcsecond they are in a perfect agreement with 
the referenced optical images both by the pulsar coordinates and 
the orientation of the PWN. For the HRC observations of 1999 we 
found a large formal systematic shift of the pulsar coordinates form 
their correct values by 22. 1 and -84.02 pixels of the CCD x and y 
coordinates, respectively. After the correction by this shift the ob- 
ject position and the PWN orientation are in agreement with those 
at the other X-ray and optical images obtained at different epochs. 

All this allows us to compare the optical and X-ray structures 
of the PWN at a subarcsecond accuracy level. To get deeper X-ray 
images of 2000 and 2006 the respective multiple data sets were 
combined making use of the merge-all v3.6 script of the CIAO 
tool. 



2.5 Wavelet filtering of the HST/WFPC2 images 

2.5.1 Wavelet filtering 

In order to study the structure of the PWN in greater detail, a filter- 
ing ha s been applied to t he image of the field around PSR B0540- 
69.3 JStarck & Murtaghl dl998l) ). This method is based on a multi- 
scale transform of the image, the isotropic undecimated wavelet 
transform (also named "a trous" wavelet transform) that produces a 
set of bands at different scales, each band having the same number 
of pixels as the image. The original image Co can be expressed as 
the sum of all these wavelet scales W: and a smoothed version of 
the image C p : 



USNO-Bl is currently incorporated inlo lire Naval Observatory Merged Aslromelric Dala set (NOMAD) which 
combine s astrometric and photometric in formation I'rom Hipparcos, Tycho-2, UCAC. Yellow-Blue6, USNO-B, and the 
2MASS rwwwarolsaravyaTtil/dai;i/rchpi\| 
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Distance (arcseconds) 
-4.6 o 4.6 



r 

C (x,y) = C p (x,y) + 2_jWj(x,y) 



(1) 



where p is the number of scales used in the decomposition. 

The filtering is obtained by applying a Hard-Thresholding to 
detect at each scale the significant coefficients : 



Wj(x,y) if Wj(x,y) is significant 
wax, >■) =„ .„ ' , .. ■ -r 

' it Wj(x, y) is not significant 



(2) 



Assuming a stationary Gaussian noise, w j(x, y) has to be com- 
pared to ko-j : 

if \Wj(x, y)\ ^ k(Tj then wj(x, y) is significant (3) 
if \Wj(x, y)\ < k<Tj then wj(x, y) is not significant 

If Wj(x,y) is small, it is not significant and could be due to noise. 
If Wj(x,y) is large, it is significant. The threshold kcrj is obtained 
by estimating the noise stan dard deviation at each sc ale cr, and 
by choosing a k value (see IStarck & MurtagrJ (2006) for more 
details). We use the FD R technique to set k in an adaptive way 
(Benjamini &Hochber3( ll995l) ). The filtered image is obtained by 
the addition of the thresholded wavelet scales Wj and the smooth 
version of the image C p . 

2.5.2 Wavelet filtering to study the optical variability of the PWN 

To study the optical variability of the PWN, we have to focus on the 
diffuse emission surrounding the pulsar. For this purpose, we have 
used a detection method to extract off the pulsar and other stars 
contained in the field. T he detection me t hod requires a multi-scale 
vision model defined in lBiiaoui & Rue ( 1995); Starck & Murtaghl 
(2006), based on the "a trous" wavelet transform described previ- 
ously. The multi-scale vision model (MVM) describes an object as 
a set of structures at different scales and a structure is a set of signif- 
icant connected wavelet coefficients at the same scale /. Assuming 
a Gaussian noise, the definition of a significant wavelet coefficient 
is given by expression {5). A prior can be introduced to change 
slightly this definition. For instance, in order to extract off the pul- 
sar and other stars, we can consider there is no interesting object 
larger than a given size. Then, we can force to 0, the wavelet co- 
efficients larger than this size. Having detected all the interesting 
objects in the fie ld (e.g., pulsar and stars) , they can be extracted by 
subtraction (see Starck & Murtagh ( 2006) for more details). The re- 
maining signal is the diffuse emission coming from the PWN. This 
is the emission we aim to study. 



3 RESULTS 

3.1 CTE effect and proper motion of PSR B0540-69.3 

Modern techniques like wavelet filtering of images helps to bring 
out faint features and filters off traces of CTE. In order to inves- 
tigate the importance of CTE effect on coordinate measurements 
or flux measurements we performed the following test on filtered 
and non-filtered images. We selected several background stars close 
to PSR B0540-69.3, including stars projected on the PWN. We 
used filtered and non-filtered F555W images for the 1995 and 2005 
epochs. We measured coordinates on both sets of images using 
I RAF imcentroid. Filtered and non-filtered images show no signifi- 
cant difference in accuracy within the image. A representative result 
obtained on filtered images is shown in Figure[2] where we plot the 
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Figure 2. Coordinate measurements of 13 stars in wavelet filtered F555W 
images for the 1995 and 2005 epochs (see text for further details.) 



distance between the pulsar and 13 stars with similar brightness as 
the filaments in the PWN, i.e., between magnitude 21 - 23, versus 
the positional uncertainty of the stars. In general, fainter stars pro- 
duce higher residuals, but the numbers are low, and there is no clear 
trend that the 2005 image give higher residuals due to the CTE ef- 
fect. We can therefore conclude that the CTE effect is unimportant 
for measurements of point sources, since the core of a star defines 
the position, which is unaffected by the CTE tail. 

Main uncertainties arise when one aligns images from differ- 
ent epochs and sets of observations. The mean accuracy for this is 
about 0.1 of a WFPC2/PC chip pixel (or a'005) (see[Z4]l. At this 
accuracy level we do not find any significant displacement of the 
pulsar position for the considered time base of 1999-2007. This is 
in agreement with the constraints on the pulsar proper motion re- 
cently reported by De Luca et al. (2007) using HST data between 
1995-2005. 

The situation becomes very different when we try to measure 
fluxes of point sources, since part of the flux is in the CTE tail. We 
have tested different aperture sizes and find that the most optimal 
size is 10 WFPC2/PC chip pixels, or ~ Of! 5 . This aperture encap- 
sulates a star of magnitude 21-23 and its CTE tail. This aperture 
was chosen for most of our photometric measurements. 



3.2 Optical photometry 

To study spectral properties of the PWN we performed detailed 
photometry of selected areas in all optical continuum filters for 
1999 and 2007. 

We selected several areas in the PWN to perform photometry 
using the optical broad band filters. The aperture radius was 10 HST 
pixels, and the apertu re was centered on th e maximum of flux of 
the "blob". Following JPe Luca et al. 120071) , we therefore centered 
on different areas for the 1999 and 2005 images. However, for the 
2007 epoch, there is no obvious dominant feature, and the aperture 
center corresponds to the assumed position of the "blob", had it 
con tinued to move with .04 times the speed of light as suggested 
by dDe Luca et al. |[2007h . An aperture of the same size was also 
placed on the apposite side of the pulsar position along the major 
symmetry axis of the PWN and at the same distance away from 
the pulsar as the "blob". We call this the "anti-blob" aperture (see 
Figure [3} . 

We also selected two 10 HST pixel apertures along the minor 
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Table 4. Results of optical spectral fits for several regions of the 0540 pulsar+PWN system. Columns 6 and 7 are the measured and dereddened fluxes with 
E(B-V)=0.20 (AV=0.62), respectively. 



Epoch 


Source 


Band 


log(Flux) 


power law 




Region Area 


name Log(v) 


observed dereddened 


a v a norm, const.* 




arcsec 2 


Hz 


pjy /jjy 


flSy 



1999 



2007 



Blob 



Blob 



14'.'3 F791W 
F547M 
F336W 



14'.' 3 



F814W 
F675W 
F555W 
F450W 
F336W 
F336W 
F336W 



14.581(15) 

14.738(8) 

14.950(13) 

14.574(18) 
14.649(12) 
14.741(21) 
14.818(20) 
14.950(13) 
14.950(13) 
14.950(13) 



1.343(18) 
1.026(20) 
0.584(82) 

1.255(20) 

1.160(22) 

0.974(26) 

0.837(26) 

-0.006 

0.188(72) 

0.510(50) 



1.496(18) 
1.275(20) 
0.989(82) 

1.403(20) 

1.357(22) 

1.224(26) 

1.152(26) 

0.398 rf 

0.593(72) e 

0.915(50/ 



1.350 



1.06C 



1.28(+3). 



1.24C 2 ) 



1999 Anti-blob 


14'.'3 F791W 


14.581(15) 


1.012(28) 


1.164(28) 




F547M 


14.738(8) 


0.707(32) 


0.955(32) 




F336W 


14.950(13) 


0.408(108) 


0.813(108) 


2007 Anti-blob 


14'.'3 F814W 


14.574(18) 


1.036(26) 


1.184(26) 




F675W 


14.649(12) 


0.976(26) 


1.173(26) 




F555W 


14.741(21) 


0.742(38) 


0.992(38) 




F450W 


14.818(20) 


0.641(36) 


0.956(36) 




F336W 


14.950(13) 


0.382(118) 


0.787(118) 



1.00Q C 



0.94(+5) 



1999 


Area 3 


14'.'3 F791W 


14.581(15) 


1.071(32) 


1.224(32) 






F547M 


14.738(8) 


0.763(34) 


1.012(34) 






F336W 


14.950(13) 


0.318(132) 


0.723(132) 


2007 


Area 3 


14'.'3 F814W 


14.574(18) 


1.086(28) 


1.234(28) 






F675W 


14.649(12) 


1.035(28) 


1.232(28) 






F555W 


14.741(21) 


0.799(36) 


1.049(36) 






F450W 


14.818(20) 


0.753(30) 


1.067(30) 






F336W 


14.950(13) 


0.276(146) 


0.681(146) 



1.30(!«) 



1.32(+«) 



0.95(«) 



1999 


Area 4 


14'.'3 F791W 


14.581(15) 


1.074(26) 


1.226(26) 


l.OOCjp 








F547M 


14.738(8) 


0.796(28) 


1.045(28) 


1.06(+3) 






F336W 


14.950(13) 


0.436(108) 


0.841(108) 






2007 


Area 4 


14'.'3 F814W 


14.574(18) 


1.098(24) 


1.246(24) 


0.86O 








F675W 


14.649(12) 


0.951(30) 


1.146(30) 








F555W 


14.741(21) 


0.746(40) 


0.996(40) 










F450W 


14.818(20) 


0.696(34) 


1.010(34) 




1.03C;!) 






F336W 


14.950(13) 


0.515(94) 


0.920(94) 







" Dots in the column show which filters were used for power law fitting. 
* The row shows which pivot frequency was used for the normalization. 
The power law was fitted without the F336W data. 



An upper limit. The aperture was placed on the estimated blob position in 2007. 
The aperture was placed on the estimated blob position in 2005. 
The aperture was placed on the estimated blob position in 1999. 



symmetry axis of the PWN, symmetrically on both sides of the 
pulsar position, and name these "Area 3" and "Area 4". Thanks to 
the high spatial resolution of the HST/WFPC2/PC chip (a' 0455 per 
pixel), we were able to measure the optical fluxes in these areas. For 
the optical background estimate for these areas we used annulus 
apertures of 6 HST pixels around 6 selected stars (see Figured- 
Photometry for all selected apertures was performed using 
IRAF/digiphot/phot. The results are presented in Tableland Fig- 
ure [4] "Area 3" and "Area 4" show slight brightening from epoch 
1999 to 2007, but with similar spectral index, a, defined as F v cc 
v~". The area "anti-blob" does not show any change within errors 
in flux and spectral index between the two epochs. 

The most interesting behaviour is displayed by the "blob". 
First of all, the flux of the "blob" became significantly lower in 



2007 compare to 1999. On other hand, the spectral index became 
shallower in 2007 (see Table [4), corresponding to a harder elec- 
tron spectrum. There are two possible hypotheses: either the "blob" 
moves with ~ 0.04c times the speed of light, as suggested by De 
Luca et al. (2007) when they compared 1995, 1999 and 2005 im- 
ages, or it fades away. We studied both these hypotheses. 

The "blob" was in its brightest phase in 1999. The aper- 
ture was centered on the maximum of the flux within the area. 
In later epochs we placed the aperture on expected positions tak- 
ing into account an assumed speed of 0.04c. For the F336W im- 
age in 2007 the expected position of the "blob" is outside the vis- 
ible PWN. Since the PWN has quite steep power law, and since 
the extinction is higher in the blue, we only have an upper limit 
on the flux of the "blob" in the F336W band (see Table [4] and 



Evolution and polarization of the pulsar wind nebula ofPSR B0540-69.3 1 




1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 




4- 

position 05 1 

position 07 ^ 



14.5 



14.7 14.8 
Log v [Hz] 



1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 




14.5 14.6 14.7 14.8 14.9 15 

Log v [Hz] 

Figure 4. Broad band optical spectra for the selected areas shown in Fig- 
ure [5] In the top panel we show the "blob" spectrum in red for 2007 and 
black for 1999. For the 2007 "blob" spectrum we have applied a 074 shift 
for the ap erture, co rresponding to a 0.04c motion of the "blob", as sug- 
gested by foe Luca et al. I l2007l) . For the spectral fit of this spectrum, we 
have omitted the F336W data (see text). The sensitivity of the F336W to 
where the aperture is placed is shown in the plot. Data in other bands are 
less sensitive (see Sect. 3.6). For the anti-"blob" we did not apply any shift 
in aperture position between 1999 and 2007. 



Figure 5. Continuum intensity (Stokes I) image at 602 nm of the PWN 
centered on PSR B0540-69.3 with overlaid linear polarization vectors. The 
size of the vector is degree of polarization in percents, while orientation of 
a vector is the position angle of linear polarization. A horizontal tick mark 
to the lower right shows 10% degree of polarization. The major axis of the 
nebula, used for Figs. [7] is marked by a solid line going from north-east 
to south-west, cutting across the pulsar and the "blob" seen here in yellow. 
The x and y axes show distances north and west of the pulsar in arcseconds. 



Figure [4j. This spectral point was not included in calculations 
of the power law for the "blob" in 2007. All power law calcu- 
lations were done using the algorithm described in Serafimovich 
et al. (2004). Fluxes in F555W and F6 75W are contaminat ed by 
emission lines of [O III] ,U4959, 5007 dMorse etal. 1 d2006h ) and 
[S II] /1/16716, 6731, respectively, which could add some system- 
atic error to the continuum flux for those data entries in Tableland 
Figure [4] 



3.3 Polarization of the PWN of PSR B0540-69.3 

Fir st optical polarizatio n meas urements for 0540 were performed 
by IChanan & Helfandl dl990h using the CTIO 4m telescope. 
Chanan & Helfand found a linear optical (V band) polarization, 
integrated over the PWN, which was P = 5.6% ± 1.0%, oriented 
at a position angle (eastward from north) 79° ± 5°. The time inte- 
grated observations with a spatial resolution of 076 per pixel did 
not allow them to detect any fine details. The first polarimetric 
phase-resolved observati ons of the PSR B0540-69.3 reported by 
(Midd leditch et al. 1119870 gave only an upper limit to the polariza- 
tion. 

A few years later, tim e integrated optical polar ization obser- 
vations were performed by Wagner & Sei fertl ( KOOOl) using the 8m 
VLT telescope. Wagner & Seifert found up to 20% of linear polar- 
ization at the rim of the diffuse nebulosity, and 5% for the pulsar. In 
this case, the spatial resolution was 072, but the seeing conditions 
(from 074 to 1") could not bring out any detailed structures. 

Polarization of the 0540 PWN has also been studied at 3.5, 
6 and 20 cm in the radio, taking into account the Faraday rotation 
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Figure 6. Polarized flux (Stokes yjQ 2 + U 2 ) as a function of optical bright- 
ness (Stokes /) over the whole PWN in 2007. Note the weak correlation 
between polarized flux and intensity. The "blob" area is the region with the 
highest intensity and smallest error bar in polarized flux. 



effect toickel etal. 1120021) . Dickel et al. obtained a total fractional 
polarization of 20% at 3.5 cm, 8% at 6 cm and 4.5% at 20 cm. 
They conclude that this large change in polarization with wave- 
length might be a result of high depolarization . The orientation is 
consis t ent with the optical value as reported by Chana n~& Helfancj 
dl990h . lDickeletal.l (2002) noted that there is a gradient in rotation 
measure, which depends on magnetic field strength and electron 
number density, from northeast to southwest. 

The latest optical polarization observations of 0540 (see de- 
scription in Section 12.3b are compatible with earlier results, but 
bring out many new interesting details due to the very high spa- 
tial resolution of HST/WFPC2. Figure[5]shows a continuum inten- 
sity image of the PWN with overlaid linear polarization vectors. 
Polarization vectors represent the electric field vector, E, which 
is perpendicular to the magnetic field. The length of the vectors 
shows the degree of polarization. The overall picture is quite com- 
plicated. As can be seen in Figure [5] the outer parts of the PWN 
have high polarization, on the order of 30% - 40%. A simi lar be- 
haviour has been seen for the Crab PWN (see lHestej d2008l) and 
references therein). 

The pulsar itself has low polarization, i.e. 5% ± 2%. This is 
consistent with earlier observations, but is in stark contrast t o re- 
cently reported results of 16% ± 4% by (Mig nani et al. |[20icl) , us- 
ing the same data set as we have used. The difference might origi- 
nate from the difficulties of separating the pulsar from the nebula, 
which has high polarization. Another source of the different results 
might be different data reduction procedures (see section |23l l, in 
particular the correction for the instrumental polarization. 

The brightest parts of the 0540 PWN in optical continuum are 
also in general the brightest in polarized light. It is highlighted by 
the correlation in Figure [6] This is perhaps not so surprising since 
the nebula emits synchrotron radiation, which can be highly polar- 
ized for ordered magnetic fields. As can be seen in Figure [5] the 
vector E is perpendicular to, or at least not aligned with the ma- 
jor symmetry axis in the bright green area. Such flow of charged 



Figure 7. Flux (Stokes /), polarized flux (Stokes ^JQ 2 + U 2 ) and polariza- 
tion angle along the major symmetry axis in 2007 (shown in Fig. [3). The 
black curve is the flux, the red curve is the polarized flux multiplied by 10, 
and the dashed line is the 3<x confidence level for the polarized flux. The 
blue horizontal straight line is the position angle of the major axis in Fig. [5] 
and the blue curve indicates the polarization angle, 6, in radians relative to 
noth (i.e, the position angle of the E vector). Note how 8 spans its full - n 
range. The pulsar stands out as a peak in the flux and polarized flux at zero 
arcseconds, and the peaks in flux and polarized flux at ~ l'/2 correspond 
roughly to the "blob" position in the 1999 and 2000 optical continuum and 
X-ray images, respectively (see Figs.[8]and|9]). 



particles, which could be a jet, may excite ejecta filaments and per- 
haps cause shock activity in the "blob" area (see section PTTt . The 
possibility of a jet al ong the major axis was briefly mentioned by 
iDe Luca etal. 1 (E007h . However, what is most surprising in Figure|6] 
is how weak the correlation is between optical continuum flux and 
polarized flux. The figure suggests that there is a general smooth 
polarized flux in the nebula, on top of which there is more spatial 
variation in the non-polarized flux. The "blob" area stands out as 
the region with the highest intensity and smallest error bar in polar- 
ized flux. Its polarized flux is slightly lower than other bright areas 
in the PWN, which may suggest less ordered magnetic field in the 
"blob region". 

To explore the polarization properties more in detail along the 
major symmetry axis, we have studied a 1" wide slice, centered on 
and along the major axis. We will later also use a similar slice for 
X-ray projections and the non-polarized optical data (see below). 
Figure [7] shows flux (black), polarized flux (red) and polarization 
angle (blue). We also show the position angle of the major axis 
(blue horizontal line). The largest maximum (around 1'.'2) to the 
right of pulsar is close to the "blob" position in 1999. We can see 
that the polarized flux jumps by a factor of ~ 3 at this position, with 
a small shift by ~ 0/.' 1 further away from the pulsar for the polarized 
flux compared to the continuum intensity. At the same time, the po- 
larization angle swings ~ 60° going from the inside of the blob to 
where the polarized flux peaks. The alignment of the E field and 
the major axis is ~ 20° at the position of peak intensity. This is 
close to what is observed for presumed shocks in extragalactic jets, 
where the E vector abruptly swings from perpendicular to the jet, 
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to parallel, at jet knots. We see a similar behavior on the opposite 
side of the pulsar at about l'.'O, although the structure is less clear. 
That the position angle of the polarization vector changes just be- 
fore the polarized flux reaches a maximum at the "blob" position, 
and perhaps close to the "antiblob" position, could point to shock 
activity (present or past) in those regions. 

We can also see two maxima in polarized flux well above 3tr 
on the same side as the "blob", but at ~ 2'.' 3 and ~ 2'.' 8 away from 
the pulsar. For the one at ~ 2'.'3, the position angle swings away 
from being parallel to the symmetry axis to being perpendicular, 
and at ~ 2'.' 8 the electric vector is parallel to that between the pulsar 
and the "blob". Note that the polarized flux is almost as bright in 
those regions as it is at the "blob" position. We see a weak, but 
similar trend, also on the opposite side of the pulsar, outside the 
"antiblob" position. 



3.4 X-ray counterpart of the "blob" 

To investigate how the X-ray emission correlates with the optical 
and polarization data, we carefully aligned the optical and X-ray 
images as described in section [2~4l and extracted intensity profiles 
along major and minor PWN symmetry axes. We chose 078 by 10" 
slit size to cover the "blob" in the best way, taking into account the 
resolution of different instruments. 

Figures[8]and[9]show an optical image in F547M taken in Oc- 
tober 1999 and an X-ray image in HRC-I taken in June 2000, with 
their respective profiles. The position of the "blob" is marked by 
vertical dotted line at the bottom panel in Fig. [8] (see also the left 
panel of Fig.[TT]for the 2006 Chandra/ ACIS data). The spatial pro- 
file clearly shows a maximum of flux from that area. As can be 
seen in the bottom panel in Fig.[9]the same optical coordinates cor- 
respond to a maximum of flux from the "blob" area also in X-rays. 
This clearly shows that the enhanced X-ray and optical radiation 
come from the same area. 

We have performed the same investigation for the second 
epoch for which we have near-coincident optical and X-ray data, 
i.e., for 2005/2006. The optical emission shows a change in mor- 
phology of the "blob" as discussed later in the paper, while the X- 
ray emission does not show any change in morphology to within the 
spatial resolution. Note that the worse spatial resolution of Chan- 
dra/ ACIS compared to Chandra/HRC. 

3.5 ACIS-S X-ray spectroscopy 

Kaaret et al. (2001) were the first to attempt a study of the X-ray 
spectral structure of the 0540 PWN. They used the observation 
carried out on August 26 1999 with the ACIS-I in a continuous- 
clocking (CC) mode. A one-dimensional image oriented approxi- 
mately along the major axis of the nebula including the pulsar was 
analyzed. Spectra extracted from the western and eastern sides of 
the pulsar, excluding a 077 region around the pulsar, showed simi- 
lar photon spectral indices of a=0.96±0. 11 and 1 . 12±0. 14, respec- 
tively, for a fixed N H column density of 4.6xl0 21 crrr 2 , assuming 
Milky Way (MW) composition of the absorbing gas. The observed 
flux was twice as bright from the western region than from the east- 
ern. Later, utilizing the ACIS-S observation obtained on November 
22-23 1999, Petre et al. (2007) measured the dependence of the 
PWN photon index on the radial distance from the pulsar. They 
extracted spectra from concentric elliptical annuli centered on the 
pulsar and aligned with the PWN major axis. A systematical soft- 
ening of the nonthermal PWN emission with radius was found, as 



is also observed in the case of the Crab PWN. However, the az- 
imuthal averaging made in the Petre et al. study does not allow to 
reveal any possible asymmetry in the photon index spatial distribu- 
tion, as suggested by the PWN morphology. A strong pileup effect 
in the ACIS-S 1999 data set precludes a relaible spectral analysis, 
owing to the conversion of two photons into a single event that has 
apparent energy equal to the sum of the photon energies, leading to 
significant spectral distortions and a loss of photons. 

The new ACIS-S data listed in Table 2 allow us to perform 
a more detailed and reliable study of the spectral structure of the 
PWN owing to a significantly smaller pileup and much longer ex- 
posures. However, the spectra from the bright PWN regions, in- 
cluding the pulsar and the blob, are still distorted by a moderate 
pileup. This leads to artificial excesses of high energy events re- 
ducing the derived photon indices. For instance, a single absorbed 
power law fit is unacceptable for the pulsar and leads to an implau- 
sible photon spectral index of only half that in previous works. 

To correct for that, the ACIS pileup model was included. To 
check its applicability we first examined the spectrum of the pulsar, 
since it is a point-like object and the pileup model was developed 
for point-like sources. We extracted data from a 3x3 pixel bin cen- 
tered on the puslar. Trailing events are seen in the ACIS images as a 
faint east-north-south-west streak, emanating from the pulsar along 
the CCD readout direction, and is caused by photons detected from 
the pulsar during frame readouts. These were placed back to the 
pulsar position using the CIAO acisreadcorr tool. The total number 
of the trailing events was estimated to be less than 1.5% of the to- 
tal count number from the pulsar region (with a count-rate of ~0.5 
counts s~'). Backgrounds were taken from a 6" radius aperture with 
the coordinates RA=05:40:36.876 and Dec=-69:22: 17.38 far away 
of the remnant. 

The spectrum was fitted in the 0.3-10 keV range with an ab- 
sorbed power-law convolved with the pileup model using standard 
XSPEC vl2.3 tools. In this test for the absorbing column density, 
we assumed MW abundances as had been done in most previous 
studies of the object. We approximated the MW abundances by the 
solar abundances of Anders & Grevesse (1989). In general, the so- 
lar metallicity according to Anders & Grevesse is 0.1 - 0.3 dex 
higher than the Milky Way abundances adopted by Wilms et al. 
(2000). However, the Anders & Grevesse abundances are widely 
used, which makes it easy to compare with previous results. As 
discussed below, our conclusions do not depend on using Anders 
& Grevesse (1989) for the Milky Way abundances. 

The resulting spectral index for the pulsar, a=0.9±0.06, is 
in excellent agreement with cf=0.92±0. 11 for the rotation-phase- 
averaged "pulsar all" spectrum of Kaaret et al. (2001) obtained 
using ACIS CC (continuous-clocking) observations, and with the 
Petre et al. (2007) value of 0.92±0.25 obtained by extracting only 
pulsar trailing events from the 1999 ACIS-S observation. In both 
these latter cases, the pileup is negligible. The fit also provides a 
reasonable PSF fraction psfrac ~ 1 treated for the pileup, and a 
plausible probablility ~0.52 that adding a photon to an event yields 
a valid X-ray event. The column density, ^ lv =(5.0±0.01)xl0 21 
cnr 2 , was compatible with that of Kaaret et al. (2001) and Petre 
et al. (2007). This ensures us that the moderate pileup in the 2006 
data-set can be reliably corrected for. The derived absorbed flux 
from the pulsar region in the 0.6-10 keV range is -1.45x10"" erg 
sm~ 2 s" 1 , which is ~ 20% lower than that of Kaaret et al. (2001). 
This may be explained by some flux contamination from the rela- 
tively bright neighboring nebular regions across the readout direc- 
tion which are collapsed together with the pulsar emission in the 
ID image in the continuous-clocking mode. It could also be caused 
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Figure 8. Top-left: ll"xll" optical image of the 0540 pul- 
sar+PWN system taken in October 1999 with the HST/WFPC2 
in the F547M band. Positions of two slices with the sizes of 
078x10" are marked by white lines. They were used for the 
extraction of the ID spatial profiles along the minor (NW-SE) 
and major (NE-SW) axes of the system shown in the top-right 
and bottom panels, respectively. Coordinate origins of the spa- 
tial axis in the profile plots coincide with the position of the 
pulsar visible as a bright point source near the PWN center. 
The position of the bright PWN "blob" seen southwest of the 
pulsar is marked by a vertical dotted line at the bottom panel. 



by a difference in assumed abundances for the Milky Way, i.e., if 
Kaaret et al. (2001) used abundances with somewhat lower metal 
content than in Anders & Grevesse (1989). 

To continue our check, we also included a pileup model for the 
emission, integrated from large extended regions of the PWN, e.g., 
from the semi-elliptic apertures encapsulating the total emission 
from the left (eastern) and right parts of the nebula along its major 
axis. For the fainter left part, where pileup is negligible, the model 
only marginally improves the fit, while for the brighter right part 
the improvements are more significant. After the pileup correction, 
the latter shows a harder spectrum than the former one, with the dif- 
ferences in a of ~ 0.2 and in the observed fluxes by a factor of ~ 2. 
That is again fully consistent with what was published by Kaaret et 
al. (2001, Table 2) based on the CC observations. This demonstrates 
that the new ACIS-S data in combination with a pileup model can 
be used also for the spectral analysis of the extended PWN, at least 
at a qualitative significance level. 

As was discussed in Serafimovich et al. (2004), MW abun- 
dances are not the best to use for 0540. To investigate that, we anal- 
ysed the pulsar spectrum using different absorbing element abun- 
dances in the MW and LMC, splitting the total column density as 
Nii=N^ w +N^ MC . The difference in absorption between MW and 
LMC gas was shown in Serafimovich et al. (2004) for the pul- 
sar/PWN system, and was then used by Park et al. (2009) in their 
study of outer regions of the SNR with the same data set. We fixed 
Ng w at the value of 0.6xl0 21 cirT 2 , while N^ MC was found from 
the fit. We obtained N^ MC =(0.8±0.02)xl0 22 cnr 2 , assuming the 
LMC abundances of Russel & Dopita (1992). The resulting spec- 
tral index was a = 0.7±0.02, which is significantly harder than in 



previous works. The derived LMC column density is about 30% - 
50% higher than the upper limit discussed by Serafimovich et al. 
(2004), and the value 0.6xl0 22 cnr 2 obtained by Park et al. (2009). 
Fixing N^ MC at the latter level produces an unacceptable fit for a 
single power-law spectrum to be typical for all young pulsars de- 
tected in X-rays. 

As discussed at length in Serafimovich et al. (2004), the LMC 
column density is unlikely to be much higher than 0.6xl0 22 citT 2 . 
To resolve this problem, we included the idea of Serafimovich et 
al. (2004) that there could be a third absorbing component, namely 
that from the ejecta of 0540. So, we fixed N^ w and N^ MC at their 
reasonable values of 0.6xl0 21 and 0.6xl0 22 cirr 2 , respectively, and 
included the third component, N^ NR , assuming that it is produced 
by the SN ejecta. It has to be dominated by heavy elements, like 
He, C, O, Si, and Fe, which are found in optical spectra of the rem- 
nant (e.g., Serafimovich et al. 2005). Since 0540 is an oxygen-rich 
SNR we started with only including hydrogen+oxygen, whereby 
we obtained a better fit (x near 1.1 per d.o. f.) with N™ R a few 
times 10" cnr 2 and O/H ~ 100 times the solar value of Anders 
& Grevesse (1989), which is 8.5 x 10~ 4 , by number. The spectral 
index remained at 0.7. 

The estimated amount of "extra" oxygen is close to that de- 
rived for SN 1987 A, which can be considered as a good model 
for an oxygen-rich event in the LMC. The model of Blinnikov et 
al. (2000) for SN 1987A was used and expanded to the size of 
0540 in Serafimovich et al. (2004). The Blinnikov et al. model 
rests on the progenitor model calculated by Nomoto & Hashimoto 
(1988) and Saio et al. (1988), which assumes that 14.7 M Q are 
ejected at the supernova explosion. The total mass on the main se- 
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Figure 9. Same as in Figure [8] but for the X-ray data taken in 
June 2000 with Chandra/HRC-I in the 0.2-10 keV band. The 
slice positions and sizes are the same as in the optical image. 
The cross in the top-left panel marks the optical position of 
the pulsar, while the vertical line at the bottom panel shows 
the position of the PWN "blob" in the optical in October 1999 
(see also the left panel of Fig.QTJfor the 2006 Chandra/ACIS 
data). As seen, despite the half year difference between the op- 
tical and X-ray observations, a likely X-ray counterpart of the 
optical "blob" can be seen at the same place as in the optical 
range. 



quence is 23 M , but several solar masses are lost in stellar winds. 
The progenitor model, which is one-dimensional, is normally re- 
ferred to as the 14E1 model (Shigeyama & Nomoto 1990). This 
model consists of well-defined shells of different elements, with 
the heavier elements toward the center. Blinnikov et al. (2000) 
used this model to calculate the broad-band filter emission from 
SN 1987A, but found that substantial mixing of the ejecta was 
needed to obtain good fits to the filter lightcurves. Figures 2 and 
3 in Blinnikov et al. (2000) show the preferred structure after mix- 
ing, and the unmixed model, respectively. Including the same ele- 
ments (H, He, C, O, Si and Fe) as in Serafimovich et al. (2004), 
the abundances relative to Anders & Grevesse (1989) for N^ NR 
were assumed to be He=5.75, C=32.5, 0=151.58, Si=65.95, and 
Fe=50.47. Other elements were neglected since they probably do 
not contribute significantly to the absorption. As a result, we ob- 
tained N™ R =(2.42±0.09)xl0 19 ctrr 2 , a=0.745±0.014, and a nor- 
malization constant of (2.30+0.25)xl0~ 3 photons cm _2 s _1 keV~' 
at the reduced x~ = 1.02 for 851 d.o.f. The harder spectrum of the 
pulsar, as compared to previous works, is likely explained by ef- 
fectively lower absorption than in case of pure MW element abun- 
dances of the absorbing matter. We note, that based on the current 
data quality and the fits statistics, we cannot justify with great con- 
fidence that the obtained splitting of the absorption in three parts 
is real, since for the pure MW and/or MW+LMC column densities 
the fits are also acceptable. However, the latter leads to an N^ MC 
value which we believe is too high (see also the discussion in Ser- 
afimovich et al. 2004). The main problem with the X-ray data is 
the pileup effect. Different pileup corrections for different element 
abundances cannot be excluded, although formally we got similar 



Table 5. Results of X-ray spectral fits for three regions of the 0540 pul- 
sar+PWN system. Photoelectric absorption along the line of sight towards 
0540 was accounted for using N% w = 0.06 X 10 22 (cm~ 2 ), N^ MC = 
0.5 X 10 22 (cm" 2 ), N S H NR = 2.42 X 10 19 (cm~ 2 ). Note, that the spectra for 
the "blob" and "anti-blob" were extracted using the same apertures as their 
optical counterparts. See text for more details. 



Region Absorbed power law model 





CI",, 


normalization 
10"* ph cm~ 2 
s -1 keV" 1 


x 2 

per dof 


pulsar 


0.74 ± 0.01 


22.97 ± 2.53 


1.01 


"blob" 


0.65 ± 0.03 


2.94 ± 0.05 


1.01 


"anti-blob" 


0.73 ± 0.04 


1.84 ±0.04 


0.93 



values of psfrac^l, and probablilityx0.5, i.e., adding a photon to 
an event yields a valid X-ray event, for all three cases of the abun- 
dances models considered. Nevertheless, including absorption from 
the supernova ejecta is physically reasonable. 

Thus, based on the results from the spectral analysis of the 
pulsar emission, we fix the MW, LMC and SNR column densities 
at the above values in the following spectral analysis of the PWN. 
To study in detail the spatial spectral variation within the PWN, we 
divided its image into rectangular spatial bins. We used 2x2 pixel 
bins for the bright central parts and 2x4 and/or 3x4 bins for the 
faint outer regions to provide enough count statistics there. Back- 
ground fluxes were taken from the same region as for the pulsar. 
The spectral data for all ACIS-S data-sets were grouped to provide 
a minimum of 15 counts per spectral bin for faint regions and of 25 



12 N. Lundqvist et al. 




Figure 10. X-ray spectral index variation map of PWN 0540. The color 
coding represents values of the spectral index a v , which varies from 0.65 
to 1.8. North is up and east to the left. The white cross is the pulsar optical 
position, and the white circle is the optical aperture centered on the "blob" 
position in 1999. Note that this region overlaps with the area which has the 
smallest spectral index, a v = 0.65 ± 0.03 (see Table |3}, in X-rays. The 
distance between the pulsar and the blob is Y.'4. 



for brighter ones. The total number of source counts varied from 
about 1000 for faint regions to about several tenths of thousands for 
the blob (with a total source count rate of ~0. 125 counts s in the 
0.5-10 keV range). The background fluxes were about 2.3 x 10 -5 
counts arcsec~ 2 s -1 , which is insignificant compared to the source 
count rates. The spectra were fitted in 0.3-10 keV range with the 
absorbed power-law. The pileup corrections were significant only 
in some relatively bright regions in the right part of the PWN, in- 
cluding the "blob". For the compact "blob", the pileup model pa- 
rameters were compatible with those of a point source (see above), 
while for other bright region, where the pileup is still significant, 
they are less constrained. The results of the fit for the blob/ant-blob 
regions are presented in the Table [5] and in Figure [KJJ we show 
the spectral index map of the PWN, compiled from the spatially 
resolved X-ray spectral analysis described above. 

As seen in Figure [TT1 the blob demonstrates the hardest emis- 
sion within the PWN. Its hardness even exceeds that of the pulsar, 
underlining the peculiarity of this structure. In general, the brighter 
the region, the harder its spectrum. This leads to the left-right asym- 
metry of the hardness of the PWN along the major axis clearly seen 
in Figure[TT] where the orientation of the NE-SW major axis slice 
is the same way as in Figures [7] [8] and [9] The hardness profile 
along the minor PWN axis is, on the other hand, almost symmetri- 
cal, and demonstrates a gradual softening of the emission with dis- 
tance from the pulsar, from a~0.9 to a > 1.6. The hardness image 
(Fig. I lOt does not reveal any hardening to extre me values for the 
assumed torus region, as seen in the Crab PWN ( Weisskop f etal. I 
2000- lMorietal. 1120041) . 



3.6 Multiwavelength spectrum of the "blob" region 

In Serafimovich et al. (2004) we studied the optical spectra for the 
1999 HST/WFPC2 epoch of several positions in what we then de- 
fined the torus and jet of the PWN. One of these positions, called 
"Area 2" in Serafimovich et al. (2004), turns out to partly over- 



lap with the "blob" discussed by De Luca et al. (2007). The re- 
sults of Serafimovich et al. (2004) hint that the optical spectrum of 
their "Area 2" in 1999 was different from other regions of the PWN 
since it could have had a somewhat steeper spectrum in the optical 
with a v = 1.58^q3q. In Table 4, the spectral slope of the "blob" has 



a v = 1.35+ 



. for the same epoch. The overlap between "Area 2" in 



Serafimovich et al. (2004) and the "blob" in Table 4 is only ~ 60%. 
The results are, however, clearly compatible even within lcr. How 
the optical emission of the "blob" varies with time, if moving, as 
well as other regions how they vary with time, is discussed in Sect. 
3.2, and listed in Tableg] 

Perhaps of more general interest is to see how the optical emis- 
sion connects to the emission at other wavelengths. In Serafimovich 
et al. (2004, their Fig. 15), a comparison between the multiwave- 
length spectrum (radio to X-rays) from the entire PWNe of the Crab 
and 0540 was made. For both objects, an extension of the X-ray 
spectrum overshoots the optical spectrum in a log(v) vs. log(F y ) di- 
agram. A further roll-off towards the radio is seen in both PWNe. 
A cubic spline fit for the entire logarithmic spectrum for the Crab 
PWN, overshoots both in the optical and radio for 0540, if fit to 
its X-ray part. We have done the same multiwavelength study here, 
but concentrating on the 0540 "blob". The optical data are again 
from 1999, but for the X-ray data, we used our reductions of the 
2006 Chandra data. These have better signal-to-noise than earlier 
data. We could have chosen optical data closer in time to the 2006 
Chandra data, but a look at Table|4]for all areas, except the "blob", 
shows that the optical flux is stable between 1999 and 2007. For 
the "blob" in Table@]we applied OM shift between 1999 and 2007, 
so the overlap between the O^SS apertures for each epoch is frac- 
tional. Nevertheless, the flux in the red HST bands for the "blob" in 
Tableware similar between 1999 and 2007 to within 5 - 10%. We 
emphasize that the X-ray spectrum of the blob has been extracted 
from the same aperture position and extent as has been used for 
the 1999 optical flux measurements. For the absorbing gas column 
density N H we used the results of Sect. 3.5. 

The multiwavelength result is shown in Fig. [12] It is clear 
that an extension of the X-ray data, considering the rather well- 
determined X-ray slope, undershoots rather than overshoots the op- 
tical data. What stands out here, like in Serafimovich et al. (2004) 
for the whole PWN, is the intrinsically steep optical spectrum com- 
pared to the X-ray spectrum. We caution that the errors calculated, 
and shown in Fig. [T2] for the optical data points, are the statistical 
errors only. If we would allow for for a modest (12-15 %) system- 
atic error of the F791W flux, one can connect all optical and X-ray 
data with a single power-law to within 2cr. The slope of that power- 
law would be a v = 0.74 ± 0.03, which is marginally different from 
the X-ray spectral slope with a v = 0.65 ± 0.03. As we noted above, 
at least part of the systematic error (5 - 10%) may come from us- 
ing the spatially slightly shifted 1999 optical data, rather than data 
closer to the epoch of the X-ray data, i.e., 2006. 

We emphasize that our finding that a single power-law may fit 
the optical/X-ray part of the spectrum for 0540 was derived only for 
the "blob", and therefore does not necessarily contradict our find- 
ings in Serafimovich et al. (2004) that spectral breaks are needed for 
the whole PWN. While we leave a discussion for the entire nebula 
for a future paper, we can make use of Figs. [U |9] and [10] to check 
whether our reductions of Chandra data make a difference com- 
pared to the X-ray spectra used in Serafimovich et al. (2004). From 
Fig.QiTJwe note that the X-ray spectral index is in the range 1.4-1.8 
in the north-west and south-east part of the PWN, and from Figs. [8] 
and [9] that the optical and X-ray fluxes are ~ 2 and ~ 3 - 5 times 
weaker there, respectively, than at the "blob" position. In Fig.ll2lwe 
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Figure 11. (left) ACIS-S X-ray intensity along the same major axis slice as in Fig[9] (right) Index spatial profile along the same slice. Coordinate origins of 
the spatial axes coincide with the pulsar position, and the vertical dotted line marks the position of the blob. Southwest is to the right of the plots. The intensity 
and index maps were smoothed with a one ACIS pixel Gaussian kernel before the sampling along the slice was made. 
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Figure 12. Multiwavelength spectrum for the "blob" region. See text for details. 
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show that a single power-law with spectral index a v = 0.74 ± 0.03 
could be enough for the "blob". Using a difference of 3 dex be- 
tween the optical and X-ray frequencies, an X-ray spectral index 
of a v = 1.6 would make an extension of the X-ray spectrum over- 
shoot the optical by more than 2 dex in flux in the north-west or 
south-east parts of the PWN. A single power-law for those regions 
is thus not possible, or actually in any part of the remnant where the 
X-ray spectral index is a v > 0.9. As a v > 0.9 is true for most of the 
PWN, a single power-law for the whole PWN will not be possible, 
and our results will in this sense not contradict the findings of Ser- 
afimovich et al. (2004). From Fig.[l0]we note that a v < 0.9 applies 
only along the major axis within ~ 2 arcseconds from the pulsar. 
The simplest explanation to this is that energy injections occur in 
this region, and that cooling of electrons take place, thus steepening 
the X-ray spectrum, when electrons move away from this region. 



4 DISCUSSION 

4.1 Morphology of the PWN 

The PWN of PSR B0540-69. 3 and SNR 0540-69 3 have been im- 
aged on several occasions. In ISerafimovich et al. I J2004h we con- 
centrated on the continuum-emitting PWN, and images showing 
the line -emitting gas of the inner ~ 4" of the SNR are nicely com- 
piled in M orse et al. I {2006). While there are some general features 
in common between continuum-emitting and line-emitting gas (like 
overall shape and concentration toward the southwest), there are 
also several important differences. 

The line-emitting gas is concentrated to the south-west, but 
with distinctions between the spatial distribution of, e.g., [O III] 
and [S II]. This can be seen in the wavelet filtered images from 1999 
shown in Fig.[T3]. The [O III] structure appears overall rounder in 
shape, whereas [S II] is more focussed to the southwestern con- 
glomerate, and shows only little emission to the northeast. A pro- 
trusion extending ~ 2" outside the main body of the emission in the 
southwestern direction is seen in both [O III] and [S II]. The struc- 
ture of this is reminiscent of the Crab chimney. What is striking is 
that the protrusion connects to the area of strongest [S II] emission 
in the remnant body, whereas the [O III] emission in the remnant 
is in general weaker in the pulsar-protrusion direction, than in the 
areas east and west of it, i.e., the areas flanking the "blob". 

The optical continuum emission from the same epoch, as ex- 
emplified by the uppermost right panel in Fig. [T3] and contours 
overlaid on the line-emission maps, is even more concentrated to 
the southwest than the [S II] emission, but the overlap between 
[S II] and the continuum is obvious. No clear indication of a protru- 
sion is seen in the wavelet filtered continuum image, which could 
be due to lower signal-to-noise. The overlap between optical con- 
tinuum and regions of enhanced [S II] emission is evidenced by 
the lower right panel of Fig.[T3] where we have made a ratio map 
of [S II] and [O III]. The highest optical continuum flux occurs 
in a conglomerate with a large [S II] / [O III] ratio. The corre- 
lation between optical continuum and the [S II] / [O III] ratio is 
indeed even more pronounced than with the [S II] emission alone. 
This could either be a physical coincidence of continuum emission 
and regions of high sulphur abundance and/or lower state of ion- 
ization, or it could signal that shock excitation is more important 
in those regions compared to regions with smaller [S II] / [O III] 
ratios. Shocked gas generally emits more [S II] than [O III] com- 
pared to photoionized gas, and the temperature of the gas is higher 
(see Williams et al. 2008). A temperature higher than the equilib- 
rium temperature in photoionization dominated plasma can also be 




Figure 13. Upper left: Wavelet filtered [O III] image of PWN 0540 obtained 
with HST/F502N in 1999 with overlaid contours from the continuum image 
in the upper right panel. Upper right: Wavelet filtered continuum image of 
PWN 0540 obtained with HST/F547M in 1999 with overlaid contours. The 
black cross shows the pulsar position. Lower left: Wavelet filtered [S II] 
image obtained with HST/F673N in 1999 with overlaid contours from the 
continuum image. Lower right: Ratio map of [S II] / [O III], with [S II] 
and [O III] from the lower left and upper left panels, respectively. Overlaid 
contours are from the optical continuum in the upper right panel. The scale 
at the bottom of the figure represents color coding for that ratio image. 



obtained in a time dependent situation. An example is the flash pho- 
toionization of the inner ring around SN 1987A where the [O III] 
temperature was ~ 5 x 10 4 K shortly after the outburst (see Fig. 1 
in Lundqvist & Fransson 1996). We explore this possibility further 
below. 

The spatial distribution of the continuum emission in the PWN 
of 0540 shows high variability in the optical. This was noticed by 
De Luca et al. (2007), and they argued that the "blob" shown in Fig. 
1 141 had moved between 1999 and 2005. The change in the spatial 
distribution of the optical continuum is clearly seen in our wavelet 
filtered images in Fig.[l5] Figure[75lalso shows how a new structure 
appears just southwest of the pulsar, as if a new "blob" would be 
generated, whereas the overall structure does not differ between the 
two epochs. In addition, Fig.[l5]shows how the structure extending 
from the pulsar (marked by a black cross) toward the "blob" con- 
tinues in the northeastern direction, although with lower intensity 
and less extent. 

The soft X-ray emission has an even more remarkable con- 
centration of emission to the "blob" position than in the optical, 
with a soft X-ray flux as bright as ~ 10% of the pulsar emission. 
The X-ray position does not seem to vary appreciably with time 
between 2000 and 2006, and the X-ray spectrum is a power-law 
which is hardest at the "blob" position, and becomes increasingly 
softer away from it (see Fig. |10|. In general, the spectrum is harder 
along the NE-SW direction than off this orientation. The optical 
continuum emission shares the general brightness distribution of 
the X-ray emission along the NE-SW direction, but it does also 
show strong variability of the spatial distribution. That we are able 
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Figure 14. The PSR B0540-69.3 and its PWN obtained with HST/F547M 
in 1992 (upper left), in 1999 (upper right), in 2005 (lower left) and 
HST/F450W in 2007 (lower right). Note that the cut levels for the F540W 
image were compensated for the slope of the nebula's SED. The white cir- 
cular apertures show the "blob" position in 1999 and respective epoch, as 
well as three other areas selected for comparison. In the 2005 and 2007 
images we have also made circular apertures for the 2005 position of the 
"blob", and the expected 2007 position if the "blob" displacement would be 
linear in time. Some stars projected on the PWN were subtracted off. 

to see this in the optical, but not in X-rays, may be an effect of 
better spatial resolution. 

4.2 Excitation of the line emitting filaments close to the 
"blob" 

The variability of the spatial distribution of the optical line emission 
is less known since observations through HST line filters have only 
been performed once, namely in 1999 (Morse et al. 2006). How- 
ever, we also have some slit spectra taken in 1996, which partly 
sample the "blob" region (N. Lundqvist et al., in preparation), and 
which do not indicate strong [S II] emission at the position of the 
1999 "blob". This, together with the fact that the optical continuum 
emission does vary, suggests that the line emission may also vary 
with time. 

The excitation of the line emitting filaments is most naturally 
explained through shock activity (Williams et al. 2008), but in Sect. 
4.1 we also highlighted the possibility of time dependent photoion- 
ization. The idea is that the optical continuum activity could also 
indicate higher X-ray activity, and that the X-rays could photoion- 
ize unshocked filaments that are quasi-stationary. As we have seen, 
the X-ray emission from the "blob" in 2006 is about a factor of 
10 weaker than the emission from the pulsar, but we only wish to 
study the photoionization close to the "blob" where X-rays from the 
"blob" dominate, so to study whether time-dependent photoioniza- 
tion is important, we have ignored the pulsar X-rays. 

To test the effect of time dependence we have made models 
where we have artificially switched on the X-ray emission from the 
"blob" to study how it would photoionize nearby gas. The model is 
similar to that used in Lundqvist & Fransson (1996) for SN 1987A, 



but with several updates (see e.g. Mattila et al. 2010). We have as- 
sumed that the "blob" is spherical with radius 072, and that the gas 
being photoionized surrounds the "blob" also in a spherical fashion. 
For the density of this gas we adopted the results of Serafimovich 
et al. (2005) who derived an electron density in the range 1 - 5 x 10 3 
cm -3 from the [S II] AA 6716-31 lines. For the abundances of the 
filaments we have assumed a composition with the relative abun- 
dances (by number) suggested by Williams et al. (2008), i.e., C: 
O: Ne: Mg: Si: S: Ar: Ca: Fe = 0.1: 1: 0.2: 0.1: 0.1: 0.1: 0.1: 0.1: 
0.1. We have also blended in H and He at various amounts as Ser- 
afimovich et al. (2005, see also Morse et al. 2006) clearly revealed 
the presence of hydrogen. The ionizing luminosity from the "blob", 
when switched on, is assumed to be 4.66 x 10 31 (£/1 keV)~ 065 erg 
s" 1 eV 1 . With H and He mixed in so that H: He: O = 1: 1: 1, 
and with a total density of atoms and ions n at = 1.5 x 10 3 cirr 3 , 
the ionization parameter (number density of ionizing photons di- 
vided by flat) is ~ 10~ 3 , i.e., only a thin rim of gas is photoionized 
around the "blob". Steady state is reached in ~ 10 years, which is 
longer than the time scale for changes in the "blob" structure, justi- 
fying time dependence. Peak temperature of the ionized gas reaches 
~ 8.5 x 10 3 K, and occurs already after 3 months. Maximum ioniza- 
tion occurs when steady state is being approached. Although sul- 
phur is mainly in S II, i.e., which we know from Figure |T3]is just 
what is observed in the "blob" region, most of the oxygen is still in 
O I. This contrasts the observations of Kirshner et al. (1989) which 
show very weak [O I] emission. One must caution that the obser- 
vations of Kirshner et al. were made more than one decade before 
the "blob" in 1999, and that the slit of Kirshner et al. was not ide- 
ally placed to probe the "blob" region. The very small ionization 
parameter and the low temperature of the ionized gas are, however, 
enough to make us rule out photoionization by X-rays from the 
"blob" structure as a viable source for the excitation of the [S II] 
emitting gas in the southwestern part of the PWN. If we disregard 
the possibility that this region is truly more abundant in sulphur 
than other parts of the nebula, the remaining explanation is what 
is discussed by Williams et al. (2008), i.e., [S II] being boosted by 
shock activity. 



4.3 Is the "blob" a moving entity? 

In Fig. [14] we also include the epochs 1992 and 2007. The signal- 
to-noise is too low in the 1992 data, which we have included for 
completeness, to draw any conclusions more than that the contin- 
uum emission is enhanced in the southwestern direction. The 2007 
image shows a much weaker "blob", and it is in fact not obvious 
that the "blob" structure is longer present in 2007. In general, the 
U-shaped structure in the 1999 image, with the pulsar at the upper 
left part of the U-structure, and the "blob" at the bottom, remains 
roughly intact between 1999 and 2005, but breaks up in 2007. A 
"blob"-structure so clearly seen in 2005 divides up into three frag- 
ments with the lower one roughly at the expected position of the 
"blob", but with the strongest emission seen northeast of the 1999 
"blob" position. In 2005 that position is actually devoid of strong 
emission. One could imagine that this new structure in 2007 could 
signal propagation of the 2005 "blob" just southwest of the pulsar, 
but we also note that the structure just northeast of the "blob" in 
2007 is clearly seen in 1999, but less brighter. A possibility is that 
the whole structure responsible for the continuum emission exists 
at all epochs, but that it contributes with different amounts at differ- 
ent epochs. With this interpretation, the northern part of the "blob" 
is just an enhancement in emission in 1999, which then fades away, 
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Figure 15. Upper: Wavelet filtered continuum image of PWN 0540 ob- 
tained with HST/F547M in 1999 with overlaid contours. Lower: Wavelet 
filtered continuum F547M in 2005 image with overlaid contours from the 
1999 epoch. Note how the "blob" structure and the structure just south- 
west of the pulsar change in morphology between these two epochs. For 
the structure close to the pulsar the strongest emission in the 2005 image 
is much wider than in the upper panel. The increased emission close to the 
pulsar could be new blob emerging. 

leaving only the southern part of the "blob" seen in 2005, and even 
more so in 2007. 

Our findings can therefore suggest that the "blob" structure 
may not be a moving entity, but rather a part of the filamentary 
structure that was excited in 1999, and then fades in optical in- 
tensity. Particularly evident is how well the continuum emission 
traces the emission structure in [S II] and how low the correla- 



tion is between continuum and [O III] (see Sect. 4.1). The high 
[S II]/[0 III]-ratio suggests shocks (cf. Sect. 4.2), this possibility 
being strengthened by the fact that the "blob" coincides with a po- 
sition in the PWN where the polarization angle changes abruptly 
and where there is strong X-ray emission. A possibility is an out- 
flow in the PWN from the pulsar which excites the filaments at the 
"blob" position in a time-varying fashion. Another possibility is a 
local release of energy at the "blob" position which then causes a 
point-like expansion of plasma. In both these cases there will be, 
after some delay, slow shocks transmitted into filaments embedded 
in the PWN in a way similar to the shocks transmitted into the blobs 
of the circumstellar ring around SN 1987A (e.g., Groningsson et al. 
2008). 

As discussed in Sect. 3.3, there are several positions along the 
major symmetry axis (NE-SW) of the PWN where the polarized 
flux peaks other than at the "blob" position, and where the polar- 
ization angle changes at these positions. This could be signs of 
past, or weak present, shock activity. A possibility is that energy 
is released along this symmetry axis at different places at different 
epochs, which means that the "blob" was just a recent such event. 
If the "blob"-like feature seen in 2005 just southwest of the pul- 
sar is a new feature moving southward, the sporadic energy release 
scenario may prove wrong. If it breaks up, it could speak in favor 
of a sporadic energy release. Indeed, we note in 2007 a peak above 
3cr at ~ ft' 6 in polarized flux, southward along the major symmetry 
axis (see Fig. 8), as well as a change in polarization angle. Look- 
ing at Fig.[l4]from 2007, no obvious blob-like counterpart is seen 
at that position. This could hint that the 2005 "blob"-like feature 
broke up, leaving only an imprint on the polarization map. 

We also note that the two strong peaks in polarized flux south- 
west of the pulsar occur in the protrusion-like feature most clearly 
seen in lower left panel of Fig. [13] which shows wavelet-filtered 
[S II] emission from 1999. The NE-SW slice in Fig.[8]catches the 
northern part of this protrusion. Further polarization studies would 
shed more light on this and other regions of the PWN. At the mo- 
ment we only have good data from one epoch, namely the one from 
2007 discussed here. 

The strong spatial variation of the inferred X-ray power-law 
slope of the emitted radiation (cf. Fig.UOIl, shows that multiwave- 
length spectra should be studied over limited spatial extents. In 
Serafimovich et al. (2004) we invoked possible spectral breaks be- 
tween the optical and X-rays when studying the multiwavelength 
spectrum of the full PWN. From Fig.[T2]here, we argue that a sin- 
gle power-law may be sufficient for the "blob". A single power-law 
could mean that the 1999 "blob" position is a site of energy injec- 
tion to relativistic electrons. Further multiwavelength, and prefer- 
ably contemporaneous, studies at other positions, away from the 
"blob", would reveal if a break in the power-law appears, and if 
cooling of the electrons occurs. We leave this for a future study. 

4.4 The supernova progenitor 

The X-ray spectral analysis presented in Sect. 3.5 gave better re- 
sults if we included "extra" absorption from oxygen-rich gas in ad- 
dition to the interstellar gas in the LMC and Milky Way. The col- 
umn density of extra oxygen we found to be needed was Nq NR as 
3.1xl0 18 cm 2 . Adopting SN 1987A, and particularly the model 
used in Blinnikov et al. (2000) and Serafimovich et al. (2004) as 
a model for 0540, N™* = 9.5xl0 17 cnT 2 for an SN 1987A age 
of t = 800 years, where the column density scales with age as 
r 2 . It seems unlikely that 0540 is younger than 800 years, so if 
2V£ is indeed as high as suggested from our X-ray analysis, then 
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more oxygen was produced in 0540 than in SN 1987A, and/or the 
oxygen-rich ejecta of 0540 expand more slowly or it could be dif- 
ferently mixed than in the model of Blinnikov et al. (2000). As we 
do not find any evidence for strong variation of the oxygen column 
density over the PWN, a higher than average column density along 
the line of sight to us is less likely. 

It is interesting to note that the model used in Blinnikov et al. 
(2000) invokes ~ 2 M Q of oxygen, whereas Williams (2010) in his 
PhD thesis argues for an upper limit of ~ 3.5 M Q of oxygen. If 
the oxygen mass is this high in 0540, and the velocity structure of 
oxygen more "compact" than in SN 1987A, Nq NR could be close to 
the value inferred from our X-ray extractions. A high oxygen mass 
would then favor a progenitor in the upper mass range for the ones 
discussed by Chevalier (2006), and we would end up with the same 
mass range of 20 - 25 M Q as favored also by Williams (2010). We 
note that Park et al. (2009) did not need to include elevated lev- 
els of oxygen for the X-ray absorption, which does not contradict 
our indications. Their analysis was made for the outer parts of the 
supernova remnant not covered by the oxygen-rich ejecta. 

4.5 Future observations 

To map out the 3D shock structure of the filaments, good reso- 
lution both in space and velocity is needed. Previous deep spec- 
troscopic observations of 0540, partly used in Serafimovich et al. 
(2004, 2005) and discussed in a forthcoming paper (N. Lundqvist 
et al., in preparation, see also above), show no clear [S II] activity 
close to the "blob" position, but these observations were made prior 
to the knowledge about the "blob", and the slits were not ideally 
placed to constrain the "blob" activity. These observations were 
also carried out at seeing > l'.'l. To map out the velocity struc- 
ture of the "blob" region in order to constrain the propagation of 
the shocks, observations with good spatial resolution and cover- 
age are needed. Modern Integral-field-units (IFUs) with a spectral 
resolution of ( < 100 km s~') are particularly useful, especially in 
good seeing and at low airmass. The good velocity resolution is 
useful to separate out slow (tens of km s~') shocks driven into the 
filaments from the faster (~ 250 km s~') shocks driven into more 
hydrogen-rich gas (cf. Williams et al. 2008). Observations in the 
optical should also be sensitive enough to catch weaker lines like 
[O III] /14363 for temperature estimates. 

The "blob" in the southwest part of the PWN, first seen in 
1999, was truly an extraordinary event. As we have argued for 
above, similar events may, however, have occurred also in the past. 
Continuous imaging monitoring of 0540 is therefore important, 
both in the continuum and in lines like [S II] and [O III], to see 
if new blobs emerge. HST can continue to play an important role 
here, especially with the upgrade to WFC3 which has no CTE 
effects. To see how polarization evolves would be particularly 
useful for a global picture of the PWN. Further down the line, both 
JWST and ALMA with their good sensitivity and spatial resolution 
will provide very important links between optical and radio studies 
of the PWN of 0540. Such studies at good angular resolution will 
also make it more easy to compare 0540 to its much closer twin, 
the Crab. 
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